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HIGH-RESOLUTION  OVER-THE-HORIZON  RADAR  USING  TIME  REVERSAL 


1.  INTRODUCTION 

Over-the-horizon  radars  (OTHRs)  [1,2]  have  the  unique  capability  of  long-range  (up  to  4,000-km) 
surveillance  over  the  horizon.  However,  OTHRs  raise  three  major  technical  issues,  as  illustrated  in  Fig.  1. 
First,  ionospheric  conditions  constantly  change,  necessitating  real-time  monitoring  and  adaptive 
beamforming.  Second,  the  resolution  cell  size  of  an  OTHR  is  typically  10  km  by  10  km,  which  is  much 
larger  than  the  size  of  targets.  Third,  there  are  many  scattering  centers  along  the  beam  path,  which  causes 
multipath  interference,  which  in  turn  causes  signal  spreading  and  inter-symbol  interference. 


Fig.  1  — Issues  of  sky-wave  HF  over-the-horizon  radar  (HF-OTHR) 


Time  reversal  has  been  successfully  demonstrated  in  acoustics  [3-7].  However,  the  implementation  of 
time  reversal  in  the  microwave  domain  has  been  impeded  by  the  lack  of  fast  enough  analog-to-digital 
converters  (ADCs).  Also,  a  skepticism  that  the  success  of  time  reversal  in  acoustics  might  not  generally 
work  with  electromagnetic  waves  has  prevailed  in  spite  of  much  success  in  optics.  Nevertheless,  time 
reversal  of  lower-frequency  electromagnetic  waves  was  demonstrated  almost  simultaneously  by  Lerosey 
et  al.  [8,9]  and  Henty  et  al.  [10].  Since  then,  time  reversal  at  electromagnetic  waves  has  been 
demonstrated  by  numerous  groups  for  various  applications,  including  communications  [11-18]. 

Phase  conjugation  (a  former  version  of  time  reversal)  has  been  considered  as  a  means  of  real-time 
compensation  of  ionospheric  distortion  without  requiring  extensive  adaptive  processing  [19].  However,  in 
addition  to  the  uncertainties  related  to  the  reciprocity  of  the  ionosphere,  the  use  of  time  reversal  for 
OTHR  has  been  impeded  by  two  other  factors:  First,  even  with  time  reversal,  the  OTHR  beam  cannot  be 
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focused  beyond  the  diffraction  limit,  which  is  typically  about  10  km  by  10  km,  and,  thus,  is  much  larger 
than  a  typical  target  size.  This  causes  significant  amount  of  clutter  (typically  50  dB  stronger  than  that 
from  a  small  target),  thus  swamping  the  target.  Second,  no  practical  means  of  imaging  method  has  been 
known  using  time  reversal.  Most  of  the  time  reversal  methods  proposed  to  date  are  mainly  intended  for 
retro-directive  beam  focusing  on  a  target  for  tracking.  Although  several  decomposition  methods  have 
been  developed  for  some  imaging  applications,  they  normally  require  a  priori  knowledge  of  the  geometry 
and  of  the  electromagnetic  properties  of  the  propagation  medium  in  addition  to  extensive  computation. 

In  this  report,  we  describe  a  new  high  resolution  method  for  OTHR  imaging  based  on  time  reversal. 
By  incorporating  the  extended  virtual  aperture  (EVA)  concept  [7,  17,  18]  enabled  by  multipath 
interference  into  a  double-pass  conjugate  imaging  scheme,  high-resolution  radar  imaging  might  be 
achieved  without  requiring  extensive  computation  or  a  priori  knowledge  of  boundary  conditions.  Also, 
proof-of-concept  experimental  results  on  several  key  features  of  time  reversal,  such  as  reciprocity, 
temporal  focusing,  spatial  focusing,  multipath-enabled  extended  virtual  aperture,  and  polarization 
diversity  time  reversal,  are  demonstrated  at  HE  or  VHP  signals. 

2.  TIME  REVERSAL  -  GENERAL 

The  following  paragraphs  describe  some  of  the  key  features  of  time  reversal  related  to  our 
experimental  results. 

2.1  Spatial  Focusing  -  Time  Reversal  Alone  Cannot  Focus  HE  OTHR 

While  time  reversal  can  be  used  to  focus  beam  on  a  target,  the  focused  beam  spot  size  is  still  limited 
by  Rayleigh  diffraction.  The  size  of  a  resolution  cell  is  given  by  5x  =  RX.  /  D  =  X.  /  [2  sin(6H)],  where  R,  X, 
D,  and  0h,  are  the  range,  wavelength  of  the  transmit  waveform,  aperture  size,  and  half  of  the  subtended 
angle  of  the  antenna  seen  at  the  target,  respectively.  For  an  HE  OTHR  with  R  =  2,000  km,  D  =  3  km,  and 
X.  =  15  m,  the  ground  resolution  5x  becomes  10  km,  which  is  much  larger  than  a  typical  target  size.  As  a 
result,  time  reversal  fails  to  focus  the  OTHR  beam  on  a  target,  making  its  usefulness  limited. 


Target 


Fig.  2  —  Time  reversal  alone  eannot  foeus  HF  OTHR  beam;  (a)  diffraetion  limit  of  eonventional  time  reversal  mirror;  (b)  the 
footprint  of  a  beam  on  target  and  sea  elutter.  TRM:  is  Time  Reversal  Mirror;  R  is  range;  D  is  aperture  size;  k  is  wavelength;  and 
Ax  is  ground  resolution. 
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2.2  Multipath-enabled  Extended  Virtual  Aperture 

In  1997,  Fink’s  group  discovered  the  effect  of  a  multipath-enabled  EVA,  which  can  have  a  significant 
impact  on  all  areas  of  imaging  [7].  Generally,  the  resolution  of  an  antenna  is  governed  by  diffraction,  and 
higher  resolution  is  obtained  with  a  larger  aperture.  However,  if  there  are  random  scatterers  along  the 
beam  path,  as  shown  in  Fig.  3,  the  aperture  size  viewed  by  the  target  can  be  increased,  and  as  a  result, 
resolution  can  be  increased  due  to  the  increased  virtual  aperture  indicated  by  the  dotted  lines.  Such  a  high 
resolution  concept  has  been  experimentally  demonstrated  by  various  groups  for  microwaves  [8-18]  and 
acoustics  [4-7]. 


Fig,  3  —  Multipath-enabled  extended  virtual  aperture.  Solid  lines  and  dotted  lines  represent  real  and 
virtual  beam  paths,  respeetively. 


Luckily,  OTHRs  operate  in  an  environment  in  which  the  various  sources  of  multipath  interference 
(such  as  ionosphere,  sea  surface,  and  ground  reflections)  are  abundant.  By  using  EVA,  resolution  of  an 
OTHR  might  be  drastically  improved  beyond  the  diffraction  limit. 

2.3  Temporal  Focusing 

In  addition  to  spatial  focusing,  time  reversal  has  the  capability  of  temporal  focusing,  which  is 
important  for  improving  resolution  along  the  range  direction  and  for  compensating  for  temporal 
distortion.  As  shown  in  Fig.  4,  a  short  pulse  (impulse)  from  a  target  (or  beacon)  is  reradiated  and 
propagated  through  a  distorting  medium.  The  waveform  experiences  distortion,  which  can  be 
mathematically  expressed  as/*  hi,  where /is  the  transmitted  waveform,  hi  is  the  impulse  response  of  the 
distorting  medium  for  the  beam  propagating  from  right  to  left,  and  *  represents  convolution.  The 
distorted  beam  is  detected,  time  reversed,  and  re-transmitted  through  the  same  medium.  Then  the  received 
signal  at  the  target  is  given  by  /  *  (/zz^/zi?),  where  hi  and  Hr  are  the  impulse  responses  for  a  signal 
propagating  leftward  and  rightward,  respectively,  and  ★  denotes  correlation.  If  the  medium  is  reciprocal, 
then  hi  =  Hr.  Moreover,  if  the  impulse  response  has  a  large  enough  bandwidth,  then  the  autocorrelation 
converges  to  a  Dirac  delta.  As  a  result,/ *  =/  In  this  way,  distortion  can  be  compensated 

without  having  to  know  the  impulse  response.  Such  a  double  correlation-based  image  reconstruction  has 
been  extensively  used  in  holographic  associative  memory  and  neural  networks  [20]. 
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Fig.  4  —  Temporal  focusing 


2.4  Time  Reversal-based  Imaging 

While  time  reversal  has  been  mainly  used  for  self-focusing  a  beam  on  a  target,  some  attempts  have 
been  made  to  obtain  images.  These  efforts  were  essentially  based  on  back-propagation  and  consisted  of 
the  following  steps.  First,  a  beam  is  transmitted  from  each  antenna,  and  the  reflected  signal  from  the 
target  plane  is  received  by  each  of  the  antenna  arrays  to  obtain  the  inter-element  impulse  responses 
(Green’s  functions).  Then  the  multistatic  matrix  K  is  formed  from  the  measured  inter-element  responses, 
and  the  time  reversal  operator  K^K  denotes  adjoint  or  complex  conjugate  transpose)  is  obtained.  The 
eigenvectors  of  the  time  reversal  matrix  are  the  singular  vectors,  and  the  eigenvalues  are  the  squares  of 
singular  values.  Each  eigenvector  is  back-propagated  through  the  medium  to  obtain  the  entire  target  field. 
Several  different  approaches  employing  the  singular  value  decomposition  (SVD)  have  been  proposed  and 
implemented.  These  include  the  DORT  algorithm  (a  French  acronym  for  the  decomposition  of  the  time 
reversal  operator)  [21],  the  MUltiple-SIgnal-Classification  (MUSIC)  algorithm  [22],  etc.  However, 
besides  extensive  computation,  these  methods  require  a  priori  knowledge  of  the  geometry  and  the 
acoustic  or  electromagnetic  properties  of  the  propagation  medium  to  calculate  the  responses  from  the 
array  elements  to  points  in  the  image  plane.  Also,  these  methods  work  mainly  with  narrow-band 
waveforms.  Further,  the  number  of  eigenvectors  is  limited  by  the  number  of  array  elements. 

More  recently,  other  time  reversal  imaging  methods,  such  as  the  direct  subtraction  of  clutter 
components  or  the  time  reversal  adaptive  interference  canceler  followed  by  time  reversal  beamformer 
(TRAIC-TRBF),  have  been  reported  [11,  12].  These  methods  also  require  prior  knowledge  of  the 
background  information  for  mathematical  back-propagation. 
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3.  DOUBLE-PASS  EVA  SCANNING  OTH  RADARS 
3.1  Overview 

Figure  5  shows  the  proposed  double-pass  time  reversal  (DP-TR)  imaging  radar  that  is  envisioned  for 
enabling  active  imaging  [21].  The  proposed  TR  radar  consists  of  a  conventional  imaging  radar  with  an 
active  beamformer  and  a  TR  mechanism  attached  to  it.  First,  the  TR  radar  steers  the  beam  along  a  certain 
direction  (Beam  1).  As  Beam  1  propagates  through  the  medium,  it  is  distorted  by  the  medium  (Beam  1’). 
The  distorted  beam  is  reflected  by  a  target  (Beam  2’)  and  returns  to  the  radar  (Beam  2).  Now  the  distorted 
waveform  associated  with  Beam  2  has  completed  the  first  round-trip  between  the  TR  radar  and  the  target. 
Next,  the  waveform  is  time  reversed,  and  a  TR  version  is  re-transmitted  (Beam  3).  Since  Beam  3  is  a  time 
reversed  version  of  Beam  2,  its  wave  front  must  closely  resemble  that  of  Beam  2.  After  propagation. 
Beam  3  becomes  Beam  3’,  which  resembles  Beam  2’,  and  is  reflected  again  by  the  target.  The  second¬ 
time  reflected  beam  from  the  target  (Beam  4’)  resembles  Beam  1’  and  its  propagated  Beam  4  must 
resemble  Beam  1,  which  is  collimated.  Therefore,  the  beam  is  coherently  focused  not  only  on  the  target 
(Beam  3’)  but  also  at  the  radar  (Beam  4),  satisfying  the  imaging  condition.  One  should  note  that  such  an 
imaging  condition  is  satisfied  regardless  of  the  distortion  caused  by  the  propagation  medium,  assuming 
that  the  medium  is  reciprocal.  The  returning  Beam  4  passes  through  the  same  beamformer  and  is  summed 
across  all  the  elements,  as  in  conventional  radars.  Then  the  steering  angle  is  incremented  to  cover  the 
entire  field  of  view. 


Distorting  Medium 


-< 

-< 

-< 


^ —  TR  Radar 
Beamformer 


3’ 


Target 


Fig.  5  —  Double-pass  scanning  time-reversal  imaging  radar 


Figure  6  shows  a  simulation  result  to  prove  the  concept  of  the  DP-TR  imaging  radar  without 
multipath  interference.  The  original  image  (shown  as  a  dotted  green  curve)  consists  of  a  single  target 
along  the  -15°  direction  and  two  closely  spaced  targets  separated  by  0.5°  around  the  0°  direction.  All 
three  targets  are  assumed  to  have  the  same  reflectance  of  0  dB.  The  reflectance  of  the  background  is 
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assumed  to  be  -20  dB  rms  (root-mean-square).  A  scanned  image  using  the  double-pass  imaging  radar 
(solid  red  curve)  clearly  shows  the  three  targets.  In  addition,  all  the  fine  details  of  the  background 
structure  are  faithfully  reconstructed.  The  only  difference  between  the  original  and  scanned  images  is  that 
the  background  clutter  noise  is  halved  in  dB.  For  example,  -20  dB  clutter  turns  out  to  be  -40  dB  after 
double-pass  time  reversal.  This  feature  increases  the  signal-to-clutter  ratio  significantly.  Further  increase 
in  resolution  is  expected  with  the  aid  of  multipath  interference  along  the  beam  path. 
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Fig.  6  —  Simulation  result  of  the  DP-TR  imaging  radar 


3.2  HF-OTHR  Beam  Focusing  due  to  Multipath 

As  explained  previously,  if  there  are  scattering  sources  along  the  beam  path,  then  the  size  of  the  point 
spread  function  can  be  reduced  significantly  by  using  the  EVA  concept.  Figure  7  shows  an  example  of 
improved  resolution  from  EVA.  Here  Zl  6^  represents  the  angle  subtended  by  the  sensor  aperture  as  seen  by 
a  target,  as  shown  in  the  figure.  Typical  sky-wave  OTHR  parameters  are  used  in  this  calculation:  a 
wavelength  of  15  m,  an  aperture  size  of  2.7  km,  a  range  of  2000  km,  and  a  clutter  reflectance  20  dB  rms 
below  that  of  the  target.  As  a  result,  the  subtending  angle  A6  becomes  approximately  0.08°  (2.7  km/2000 
km  =  1.35  X  10'^  radians  =  0.08°).  One  should  differentiate  this  subtending  angle  from  angular 
resolution,  which  is  approximately  half  a  degree  for  the  illustrated  OTHR  case.  With  a  conventional  sky- 
wave  OTHR  without  EVA  (Fig.  7(a)),  the  resolution  cell  size  is  large,  approximately  10  km  by  10  km. 
Therefore,  the  signals  from  small  targets  that  are  approximately  100  m  long  are  swamped  by  the  clutter 
and  are  thus  not  detectable  without  resorting  to  other  means,  such  as  Doppler  filtering. 
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Fig.  7  —  Resolution  of  sky-wave  OTHR:  (a)  without  EVA;  and  (b)  with  EVA  and  a  subtending  angle  of  +  5° 


By  using  EVA  with  an  increased  subtending  angle  of  ±5°  (Fig.  7(b)),  the  resolution  cell  size  along  the 
cross-range  direction  becomes  as  small  as  80  m,  more  than  a  100  times  improvement  in  resolution.  As  a 
result,  most  of  the  beam  can  be  focused  on  a  target,  so  that  the  target  signal  becomes  much  stronger  and 
sharper  compared  to  the  case  without  EVA.  It  should  be  noted  that  the  resolution  is  determined  by  the 
subtending  angle,  regardless  of  the  locations  of  the  scatterers  along  the  range  direction.  Also,  the 
distribution  of  the  scatterers  can  be  sparse,  but  the  separation  between  adjacent  scatterers  must  be  random 
to  suppress  the  grating  lobes. 

Figure  8  shows  the  resolution  as  a  function  of  half  of  the  subtending  angle,  which  is  given  by  5x  =  / 

[2sin(A6^/  2)],  where  5x,  X,  and  A  6^  represent  resolution,  wavelength,  and  subtending  angle,  respectively. 
The  curve  starts  from  oo  at  0°,  rapidly  decreases  at  small  angles,  and  slowly  converges  to  half  the 
wavelength  as  the  angle  is  increased  to  90°.  As  explained  previously,  current  HF  OTHRs  operate  at  half 
of  the  subtending  angle  of  0.04°,  which  gives  a  10-km  resolution.  If  the  subtending  angle  A/9is  increased 
to  10°  (±  5°),  for  example,  the  resolution  cell  size  can  be  as  small  as  80  m  as  shown  with  the  pair  of 
dotted  lines. 

In  order  to  achieve  super-resolution  via  EVA,  a  large  number  of  randomly  distributed  scatterers  are 
required  along  the  beam  path  to  suppress  sidelobes  efficiently  along  unwanted  directions.  Also,  the 
system  is  expected  to  be  sensitive  to  external  noise  perturbations. 

Along  the  range  direction,  resolution  is  given  by  Ar  =  c  /  (2Af),  where  c  and  Af  are  the  speed  of  light 
and  the  frequency  bandwidth,  respectively.  Currently,  a  15  kHz  bandwidth  is  used  to  yield  a  10-km  range 
resolution.  If  the  bandwidth  can  be  increased  to  1.5  MHz,  range  resolution  can  be  improved  to  100  m.  In 
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this  way,  the  resolution  can  be  improved  along  both  the  range  and  azimuthal  directions,  permitting  a  time 
reversed  HF  OTHR  to  track  small  targets  with  a  high  signal-to-clutter  ratio. 


Fig.  8  — ^Resolution  of  the  HF  OTHR  as  a  funetion  of  half  of  the  subtending  angle  (A0  /  2) 


3.3  Conjugate  Imaging 

In  a  typical  sky-wave  HF  OTHR  case,  the  clutter  noise  is  50  dB  stronger  than  the  target  return  signal, 
as  can  be  seen  from  Table  1,  assuming  that  the  target  size  is  10  m  x  10  m,  target  reflectivity  is  1,  and  the 
clutter  reflectivity  is  0.1.  Consequently,  at  each  antenna  element,  the  signal  from  the  background  is  50  dB 
stronger  than  that  from  a  small  target,  so  that  it  is  not  detectable  without  Doppler  processing.  As  a  result, 
it  appears  that  the  extraction  of  the  target  signal  from  the  huge  background  is  extremely  difficult,  no 
matter  what  time  reversal  might  do. 


Table  1  — RCS  of  Target  and  Clutter  for  a  Typical  Sky-Wave  HF  OTHR 


Target 

Clutter 

Area 

lOmx  10  m 

10  km  X  10  km 

Reflectivity 

1 

0.1 

RCS 

10^  (20  dB-sm) 

10^  m"  (70  dB-sm) 

Generally,  the  output  image  i{x)  from  an  imaging  system  is  given  by  5(x)*/z(x),  where  h(x)  is  the  point 
spread  function  (PSF)  of  the  imaging  system  and  *  denotes  convolution  as  before.  As  a  result,  the  original 
image  is  smeared  by  the  size  of  the  PSF.  Here,  one  should  note  that  the  beam  from  outside  the  PSF  does 
not  contribute  to  its  corresponding  image  point.  This  observation  can  be  more  fully  appreciated  by 
analogy  to  a  camera  system,  which  forms  a  sharp  image  of  a  small  point  surrounded  by  large  bright 
background.  Although  the  total  sum  of  light  from  the  background  is  much  stronger  than  that  from  a  small 
point,  one  can  still  clearly  see  the  spot  at  the  image  plane,  without  being  swamped  by  the  strong 
surrounding  background. 
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In  the  proposed  system,  such  an  imaging  condition  is  automatically  satisfied  regardless  of  the 
intermediate  distorting  medium.  Also,  due  to  the  EVA  employed  in  the  system,  the  PSF  can  be  made  very 
small  as  shown  in  Fig.  7(b).  In  addition,  the  EVA  allows  a  sharp  beam  focusing  selectively  on  a  target. 
This  feature  of  high-resolution  beam  focusing  on  both  target  and  imaging  radar  allows  high  a  signal-to- 
clutter  ratio  and,  thus,  small  target  detection,  even  when  the  clutter  power  is  much  greater  than  the  target 
return  power. 

3.4  Signal-to-clutter  Ratio  with  Scanning 

By  incorporating  multipath-enabled  super-resolution,  which  greatly  reduces  the  size  of  the  point 
spread  function,  the  signal-to-clutter  or  signal-to-noise  ratio  can  be  significantly  improved,  as  explained 
previously.  Furthermore,  in  our  proposed  scheme,  the  beam  is  scanned  sequentially  along  both  directions 
as  with  conventional  phased  array  radar  imaging.  Therefore,  the  instantaneous  field  of  regard  of  the 
system  is  very  small,  thus  relaxing  the  difficulty  of  discriminating  the  signal  from  the  background  noise. 


Fig.  9  —  Various  issues  of  TR  HF  OTHR 


4.  TECHNICAL  ISSUES 

The  proposed  time-reversed  HF  over-the-horizon  radar  (TR  HF  OTHR)  raises  a  lot  of  technical  issues 
(Fig.  9),  which  are  described  in  this  section. 

4.1  Reciprocity 

Time  reversal  is  based  on  the  assumption  that  the  propagation  medium  is  reciprocal.  For  ionospheric 
propagation,  two  factors  can  invalidate  the  reciprocity  theorem.  First,  time- varying  ionospheric 
inhomogeneities  can  cause  the  Green’s  function  between  two  points  to  change  during  the  transit  time  of 
the  beam.  Generally,  the  drift  of  the  ionosphere  is  slow,  with  typical  drift  velocities  on  the  order  of  100 
m/s,  compared  with  the  round  trip  time  (10  ms).  Second,  the  forward  and  backward  Green’s  functions  can 
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be  different,  due  to  the  asymmetry  of  the  dielectric  tensor  or  to  some  nonlinearity.  Faraday  rotation 
induced  by  the  geomagnetic  field  is  non-reciprocal  in  this  sense. 

Time  reversal  works  if  the  all  of  the  components  of  a  scattered  radar  signal  propagate  along  a  single 
reversible  path.  When  an  HF  wave  enters  the  ionosphere,  it  splits  up  into  its  characteristic  polarizations. 
Each  polarization  propagates  at  a  distinct  phase  and  group  velocity.  For  quasi-longitudinal  propagation 
with  the  wave  normal  within  80  degrees  of  the  ambient  magnetic  field  B,  these  polarizations  are  left-hand 
and  right-hand  circular. 

A  sample  ray-trace  (Fig.  10(a))  shows  that  the  two  characteristic  modes  in  the  ionosphere  take 
different  paths  between  two  points  on  the  ground.  The  O-mode  and  X-mode  polarizations  are  LHCP  and 
RHCP,  respectively.  It  they  are  excited  with  equal  amplitude  at  the  transmitter  in  the  form  of  a  linearly 
polarized  wave,  the  electric  field  vector  at  the  end  point  will  have  rotated.  If  only  one  mode  is  excited  at 
the  transmitter,  it  can  retain  the  same  polarization  along  the  propagation  path. 

Usually,  the  waves  retain  their  mode-identity  as  they  propagate  in  the  ionosphere  with  weak  coupling 
(Fig.  10(b)).  In  an  irregular  ionosphere,  each  characteristic  polarization  may  couple  into  the  other  mode 
near  regions  where  propagation  is  transverse  to  B  (Fig.  10(c)).  Fortunately,  this  mode  coupling  is  rare. 
When  the  radar  signal  scatters  from  a  ground  target,  a  circular  polarization  of  one  mode  can  be  converted 
into  the  circular  polarization  of  both  modes.  The  reflected  ray  paths  are  subjected  to  propagation  effects 
of  the  transmitted  waves  so  that  a  single  HF  radar  transmission  can  have  many  propagation  paths.  This 
effect  is  minimized  by  launching  the  wave  with  a  characteristic  polarization  of  the  ionosphere  and 
receiving  the  scattered  wave  with  an  antenna  sensitive  only  to  that  particular  polarization.  The 
retransmission  after  time  reversal  then  launches  the  wave  that  was  confined  to  one  propagation  mode 
(Fig.  10(d)). 

Also,  for  broadband  signals,  the  transfer  function  of  an  antenna  is  different  between  transmit  and 
receive  modes.  Kanda  [24]  has  shown  that  the  transmitting  transient  response  of  an  antenna  is 
proportional  to  the  time  derivative  of  the  receiving  transient  response  of  the  same  antenna,  thereby 
causing  an  extra  ’yco"  factor  in  the  antenna  reciprocity  relationship.  This  additional  factor  can  have  a 
dramatically  different  effect  on  differentiation  and  integration,  depending  upon  the  usage  of  an  antenna 
[25,  26]. 

In  the  case  of  HF  OTHR,  the  failure  of  reciprocity  for  ionospheric  paths  has  been  reported  previously 
in  the  literature  [27-30].  However,  the  results  of  lull  and  Budden  [31-33]  suggest  that  the  effects  of 
ionospheric  non-reciprocity  can  be  significantly  reduced  by  carefully  choosing  the  polarizations  of 
transmit  and  receive  antennas.  One  means  for  effecting  this  reduction  is  to  use  polarization-diverse 
antennas.  When  non-reciprocity  is  still  an  issue,  our  double-pass  time  reversal  scheme  has  the  potential  to 
circumvent  the  problems  related  to  non-reciprocity  at  the  cost  of  extra  travel  time  for  the  double  pass. 
Furthermore,  the  temporal  non-reciprocity  also  might  be  compensated  based  on  a  priori  knowledge  of  the 
behavior  of  ionospheric  movement. 
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13.  MHz  RayPaths,  foF2  =  6.5  MHz 
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Fig.  10  —  (a)  0-Mode  (red)  and  X-Mode  (blue)  rays  between  transmitter  and  point  on  the  ground.  The  differenee  in  phase 
veloeity  between  the  two  modes  produees  Faraday  Rotation,  (b)  Weak  eoupling  in  the  ionosphere  for  waves  erossing  a 
transverse  region  in  the  plasma,  (e)  Strong  eoupling  in  the  ionosphere  for  waves  erossing  a  transverse  region  in  the  plasma,  (d) 
Optimum  eontrol  of  polarization  effeets  by  seleetive  generation  and  reeeption  of  HF  radar  signals. 


4.2  Multipath 

The  proposed  system  strongly  relies  on  the  existence  of  multipath  in  the  propagation  medium.  The 
ionospheric  environment  provides  various  sources  of  multipath  from  anisotropy,  distortion,  dispersion, 
multiple  hops,  multiple  layers,  multiple  modes,  etc.  These  multipath  effects  contribute  to  beam  focusing 
in  both  time  and  space.  In  fact,  multipath  spreads  the  temporal  signals,  and  the  temporal  impulse  response 
becomes  random,  acquiring  a  large  bandwidth  so  that  its  auto-correlation  converges  to  a  Dirac  delta,  as 
shown  in  Fig.  4.  Also,  multipath  contributes  to  high-resolution  spatial  focusing,  due  to  the  EVA  concept 
described  previously.  Therefore,  the  existence  of  a  large  number  of  randomly  distributed  multiple 
scattering  sources  is  a  pre-requisite  for  a  successful  spatio-temporal  focusing  and  EVA. 

One  might  be  concerned  that  the  ionosphere  itself  may  not  provide  sufficient  multipath  to  allow 
super-resolution.  However,  the  scattering  sources  do  not  need  to  exist  inside  the  ionosphere.  They  could 
be  near  the  target  -  for  example,  the  HF  beam  might  propagate  without  any  beam  deviation  until  it  arrives 
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near  the  target.  The  scattering  sources  near  the  target,  such  as  mountains,  buildings,  or  even  the  clutter 
itself,  also  may  contribute  as  sources  of  multipath,  as  long  as  they  help  increase  the  subtending  angle. 

4.3  Coherence,  Bandwidth,  Phase  Noise,  and  Synchronization 

Since  the  time  reversal  process  is  based  on  the  coherent  summation  of  scattered  beams  along  different 
paths,  it  is  essential  to  maintain  coherence  both  in  space  and  time,  even  for  the  longest  delayed  multipath 
signals.  Accomplishing  this  coherence  requires  a  very  stable  signal  source  with  extremely  low  phase 
noise  over  the  coherent  integration  time.  To  ensure  that  the  scattered  beams  along  different  path  lengths 
are  recombined  in  a  coherent  manner,  the  coherence  length  of  the  signal  source  must  be  sufficiently  long 
compared  to  the  variation  of  the  path  lengths. 


It  has  been  shown  that  the  peak-to-noise  ratio  (in  amplitude)  varies  roughly  as  (MAv  /5v  where 
Av  is  the  bandwidth  of  the  initial  emitted  pulse,  M  is  the  number  of  antenna  elements,  and  5v  defines  the 
correlation  frequency  of  the  reverberated  field,  which  is  roughly  the  inverse  of  the  spread  time.  The 
spread  time  is  defined  as  the  RMS  delay  spread  of  the  signal  [9]. 

The  frequencies  available  to  an  HF  OTHR  are  limited  by  regulations  of  the  Federal  Communications 
Commission  and  the  National  Telecommunications  and  Information  Administration,  as  well  as  by  the 
bandwidth  constraints  imposed  by  the  ionosphere.  For  example,  polarization  changes  associated  with 
Faraday  rotation  produce  a  variation  of  the  received  signal  amplitude  and  thus  may  impose  a  bandwidth 
limitation  on  time  reversal. 

Finally,  the  distribution  of  a  suitable  timing  source  for  synchronization  and  coherent  operation  is  one 
of  the  most  challenging  subjects  in  the  implementation  of  such  a  coherent  radar  array  [34]. 

4.4  Mutual  Coupling 

Mutual  coupling  can  cause  serious  noise  problems  in  an  antenna  array  and  with  targets.  However,  it  is 
anticipated  that  the  adaptive  nature  of  time  reversal  might  automatically  account  for  these  undesired 
effects. 

4.5  Compatibility  with  Doppler  Filtering 

Doppler  filtering  in  HF  OTHRs  plays  a  significant  role  in  removing  clutter  noise.  Therefore,  it  is 
highly  desirable  to  maintain  this  capability  even  with  time  reversal.  In  this  regard,  one  might  wonder 
whether  the  time  reversal  process  is  compatible  with  Doppler  filtering.  The  first  issue  is  whether  time 
reversal  could  negate  Doppler  frequencies.  Since  the  time  reversal  process  simply  reverses  the  timing 
signal,  it  does  not  affect  the  signal  spectrum  (frequencies).  The  second  issue  is  whether  Doppler  filtering 
could  wash  out  the  scattered  beam  from  multipath  scatterers.  To  implement  an  EVA,  multipath  scattering 
is  required,  as  explained  previously.  These  scattering  sources  are  mostly  stationary  or  slowly  moving 
objects,  and  their  return  signals  will  be  separated  out  by  Doppler  processing.  However,  one  should  note 
that  the  portion  of  the  scattered  beam  which  is  re-scattered  from  the  target  is  Doppler  shifted.  It  is  that 
portion  of  the  beam  that  contributes  to  the  EVA.  Since  its  frequencies  are  Doppler-shifted  by  the  fast- 
moving  targets,  it  will  not  be  separated  out  by  Doppler  processing. 

5.  EXPERIMENTAL  RESULTS 

Since  there  are  no  previously  reported  results  on  time  reversal  at  HF  or  VHF  to  the  best  of  our 
knowledge,  there  are  several  technical  issues  such  as  reciprocity  (especially  in  the  ionosphere),  spatio- 
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temporal  focusing,  extended  virtual  aperture,  etc.  Also,  the  proposed  conjugate  imaging  concept  must  be 
verified.  Here  we  demonstrate  several  experimental  results  to  address  these  issues. 

5.1  Reciprocity  at  HF  (Non-ionospheric) 

Figure  11  shows  a  schematic  diagram  of  the  setup  to  characterize  reciprocity  and  temporal  focusing 
of  HF  ground  wave  at  a  short  distance  (approximately  10  m).  Each  test  consists  of  two  cycles:  in  the  first 
cycle,  an  impulse  signal  is  transmitted  and  received,  and  in  the  second  cycle,  the  received  signal  is  time 
reversed,  retransmitted  and  received. 


Arbitrary  Waveform 
Generator 


Filter  /  Amplifier 


■ 


¥ 


4 


<mppi> 


Fig.  11  —  Experimental  setup  for  reeiproeity  and  temporal  foeusing  at  a  short  distanee  with  20  MHz  HF  ground  wave 


All  the  measurements  are  performed  in  the  time-domain  in  an  all-digital  manner,  unlike  most  of  other 
time  reversal  experiments  done  in  frequency  domain  with  the  aid  of  analog  components  such  as  vector 
modulators  or  down  converters  [8-15].  Such  an  all-digital  time-domain  implementation  permits  instant 
and  flexible  operation  of  an  OTHR  system.  Since  the  bandwidths  of  HF  signals  are  relatively  narrow, 
analog-to-digital  converters  and  digital  oscilloscopes  are  available  at  that  speed. 

The  initial  impulse  (represented  by  blue  solid  lines)  is  a  20  MHz  mono-cycle  pulse  (pulsewidth:  100 
ns)  generated  by  an  arbitrary  waveform  generator  (Tektronix,  model  AWG5014).  The  signal  is  then 
connected  to  a  monopole  antenna  (Hy-Gain,  model  AV-18VS)  after  passing  through  an  amplifier  and  a 
harmonic  filter.  The  transmitted  signal  is  then  detected  by  a  receiving  antenna,  which  is  a  Nearly 
Vertically  Incident  Skywave  (NVIS)  antenna  (B&W,  model  BWD-65).  The  NVIS  antenna  is  a  folded 
dipole  antenna  that  can  cover  a  broad  frequency  range  of  3  to  30  MHz.  This  NVIS  antenna  was  chosen  for 
future  skywave  experiments  at  a  short  distance. 
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The  received  signal  is  then  digitized  by  a  digital  oscilloscope  and  time  reversed  by  a  computer.  The 
time-reversed  waveform  is  loaded  on  the  arbitrary  waveform  generator  and  is  sent  to  the  monopole 
antenna  again  (represented  by  red  dotted  lines).  The  signal  received  by  the  NVIS  antenna  the  second  time 
is  detected  and  digitized  by  the  oscilloscope. 

Figure  12  illustrates  some  reciprocity  data  at  HF.  At  first,  the  monopole  antenna  is  used  as  a 
transmitting  antenna,  and  the  NVIS  antenna  is  used  as  a  receiving  antenna.  The  detected  timing  signal  is 
represented  by  a  blue  line.  Then  the  antennas  are  interchanged  so  that  the  NVIS  antenna  is  used  for 
transmission  and  the  monopole  antenna  is  used  for  reception.  The  detected  signal  for  this  configuration  is 
represented  by  the  red  line.  The  green  curve  shows  the  input  pulse.  The  data  shown  here  are  a  randomly 
chosen  sample  of  real-time  raw  data  without  any  signal  averaging  or  normalization.  Yet,  one  can  see  that 
both  curves  coincide  reasonably  closely,  at  least  in  phase,  even  though  the  antennas  are  different  in  terms 
of  gain  and  polarization.  This  result  clarifies  that  there  are  no  significant  non-reciprocity  effects  at  HF 
other  than  that  originated  from  the  ionosphere. 
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Fig.  12  -  Reciprocity  between  two  different  antennas  at  a  HF  ground  wave 


5.2  Temporal  Focusing  of  HF  Signals 

Figure  13  shows  experimental  data  demonstrating  temporal  focusing  of  an  HF  ground  wave.  Figure 
13(a)  is  the  returned  signal  in  the  first  cycle,  as  shown  in  Fig.  11.  Figure  13(b)  is  an  experimental 
correlation  result  obtained  from  the  re-transmitted  time  reversed  signal  in  the  second  cycle.  A  good  trend 
of  temporal  focusing  is  exhibited  even  with  a  small  number  of  random  scatterers  between  the  transmit  and 
receive  antennas.  Figure  13(c)  is  the  calculated  auto-correlation  of  the  signal  shown  in  Fig.  13(a).  One  can 
see  good  agreement  between  the  calculated  and  experimentally  obtained  correlation  results.  This  is 
expected  because  ground  wave  propagation  is  predominantly  reciprocal  in  nature,  unlike  sky  wave 
propagation,  which  is  susceptible  to  the  temporal  variation  of  ionospheric  irregularities  or  non-reciprocal 
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Faraday  rotation  due  to  the  magneto-ionic  effects.  Yet,  the  results  again  verify  that  reciprocity  exists  at 
HF  between  the  two  antennas. 
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Fig.  13  -  Temporal  focusing  of  an  HF  ground  wave:  (a)  first-time  returned  signal;  (b)  experimental  temporal 
focusing  after  time  reversal,  retransmission,  and  reception;  and  (c)  computed  autocorrelation  of  the  signal  in  (a).  The 
experimental  results  in  (b)  are  in  good  agreement  with  the  simulation  result  in  (c). 


In  order  to  obtain  a  desired  correlation  similar  to  a  delta  function,  which  is  the  impulse  response  of  a 
time  reversed  system,  the  one-way  impulse  response  signal  must  satisfy  the  following  requirements:  it 
must  be  reciprocal  in  both  directions,  random  with  large  bandwidth,  and  bipolar  with  a  zero  average.  This 

focusing  feature  can  be  further  enhanced  by  an  increase  in  the  number  of  antennas  in  proportion  to  Vy  , 
where  N  is  the  number  of  antennas. 

5.3  Spatial  Focusing  and  Extended  Virtual  Aperture  Experiment  at  250  MHz 

To  demonstrate  the  spatial  focusing  feature  of  time  reversal,  we  constructed  an  indoor  experimental 
setup  scaled  down  to  250  MHz  (wavelength  =  4  ft)  in  a  laboratory.  Figure  14  shows  a  schematic  diagram 
and  photos  of  the  setup. 
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Fig.  14  -  An  experimental  setup  to  test  spatial  foeusing  and  extended  virtual  aperture  at  250  MHz.  A  sehematie  diagram  (top) 
and  photos  (bottom)  of  eaeh  eomponent.  (a)  four  antennas  forming  a  time  reversal  mirror  (TRM);  (b)  eontrol  instrument;  (e) 
beaeon  /sean  antenna;  (d)  and  (e)  are  multipath  generators. 


Home-made  dipole  antennas  tuned  to  250  MHz  are  used  in  this  experiment.  Four  antennas  in  a  row 
(Fig.  14(a))  at  two-foot  (half  the  wavelength)  spacings  are  used  to  form  a  time  reversal  mirror.  Figure 
14(b)  shows  the  instruments  used  in  this  experiment.  A  digital  oscilloscope  and  an  arbitrary  waveform 
generator  (same  as  described  previously),  power  sensor  (Agilent,  U2000A,  USB  Average  Power  Sensor) 
to  scan  the  beam  pattern,  and  a  rubidium  frequency  standard  (SRS,  model  FS725)  that  was  used  as  an 
external  reference  to  both  the  oscilloscope  and  the  AWG.  All  the  signal  flows  are  controlled  with  a 
MatLab  code.  One  antenna  on  a  cart  (Fig.  14(c))  is  used  as  a  beacon  and  a  scanning  antenna. 

The  system  is  very  sensitive  to  any  electromagnetic  interference,  human  interruption,  cable 
arrangements,  or  ground  connection.  In  order  to  obtain  reproducible  results,  special  care  was  taken  to 
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avoid  these  sources  of  interference.  Also,  a  couple  of  ferrite  beads  are  added  on  each  coaxial  cable. 
Further,  several  RF  foam  absorbers  were  used  to  block  some  unwanted  reflections.  The  absorber  was  a 
rubberized  foam  material  impregnated  with  carbon  and  iron.  The  material  was  borrowed  from  another 
laboratory  and  was  not  optimized  in  any  way.  We  used  four  24-in.  by  24-in.  panels  that  were  ridged  and 
four  24-in.  by  24-in.  panels  that  were  square-cross-section  cones. 

In  order  to  generate  multipath  interference,  several  metallic  carts  covered  with  aluminum  foil  (Fig. 
14(d))  or  several  wires  hanging  from  a  clothesline  (Fig.  14(e))  were  used.  We  stretched  a  cord  across  the 
width  of  the  laboratory,  centered  between  the  scan-antenna  and  the  array.  The  height  of  the  cord  was  the 
same  as  the  tops  of  the  vertically  oriented  dipoles.  An  array  of  insulated  wires  23  in.  long  (approximately 
half  the  wavelength)  was  taped  to  the  cord.  The  wires  were  placed  along  the  cord  over  a  distance  of 
approximately  15  ft.  We  had  previously  demonstrated  that  the  hanging  wires  had  a  significant  effect  on 
the  RF  environment. 

In  this  spatial  focusing  experiment,  a  two-step  process  was  used,  starting  in  several  places  on  the  x- 
axis. 


•  Step  1:  Transmit  a  250  MHz  CW  signal  from  the  cart-mounted  beacon  antenna  located  at  a 
certain  point  along  the  x-direction.  Receive  the  signal  with  the  four-antenna  array,  record, 
and  perform  time  reversal  with  a  computer. 

•  Step  2:  Transmit  the  time-reversed  signal  from  the  four-antenna  array  in  a  synchronous 
manner.  Scan  the  beam  pattern  with  the  cart-mounted  antenna  along  the  x  direction. 

Figure  15  shows  an  experimental  result  demonstrating  spatial  focusing  and  extended  virtual  aperture 
due  to  multipath  interference.  Figure  15(a)  shows  the  result  without  multipath,  and  Fig.  15(b)  shows  the 
result  with  multipath.  The  arrows  on  the  top  represent  the  beacon  locations  for  the  corresponding  curves. 
One  can  see  that  the  beam  is  focused  at  the  corresponding  beacon  location.  The  dotted  pink  curve  shows 
the  theoretical  beam  pattern  when  the  beacon  is  located  at  the  origin.  As  is  shown  in  Fig.  15(b),  the 
beamwidths  with  multipath  are  significantly  narrower  than  the  theoretical  diffraction  limit  (shown  with  a 
dotted  violet  curve  in  Fig.  15(a)).  As  the  number  of  multiple  scatterers  increases,  the  beam  peaks  will 
become  sharper  and  narrower  permitting  super-resolution. 
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Fig.  15  -  Experimental  spatial  foeusing  and  extended  virtual  aperture  at  250  MHz:  (a)  without  multipath  (free  spaee);  and  (b) 
with  multipath  interferenee.  Dotted  violet  eurve  shows  theoretieal  beam  pattern  when  the  beaeon  is  loeated  at  the  origin.  The 
arrows  on  the  top  represent  the  beaeon  loeations  for  the  eorresponding  eurves.  One  ean  see  that  the  beam  is  foeused  at  the 
eorresponding  beaeon  loeation.  Furthermore,  the  beamwidths  with  multipath  are  signifieantly  narrower  than  the  theoretieal 
diffraetion  limit. 


5.4  Polarization  Diversity  Time  Reversal 

As  described  previously,  spatial  ionospheric  non-reciprocity  effects  can  be  significantly  compensated 
using  polarization  diversity  antennas.  Although  a  good  experimental  ionospheric  model  to  simulate  the 
anisotropic  non-reciprocity  is  not  available  at  the  moment,  we  want  to  make  sure  that  time  reversal  works 
properly  with  two  orthogonal  polarizations  -  horizontal  and  vertical.  When  there  is  no  multipath 
interference  (shown  in  Fig.  16(a)),  the  horizontal  polarization  component  shows  a  weak  gain  (about  20% 
of  the  vertical  polarization  component).  This  is  attributed  to  the  geometry  of  the  experimental  setup.  With 
multipath  interference,  as  shown  in  Fig.  16(b),  the  horizontal  components  are  significantly  increased  due 
to  the  polarization  scrambling  caused  by  the  multipath  interference.  All  these  curves  exhibit  the  proper 
spatial  focusing  feature  at  the  beacon  location. 

Currently,  a  HF  skywave  test  site  is  being  built  between  the  Chesapeake  Bay  Detachment  (CBD)  and 
Tilghman  Island  in  Maryland. 

6.  CONCLUSIONS 

We  have  proposed  a  new  OTHR  concept  using  time  reversal  that  can  address  the  three  major  issues  of 
current  OTHR  systems.  Also,  we  have  proposed  a  time-reversal-based  imaging  method  that  does  not 
require  detailed  information  on  the  boundary  conditions.  The  proposed  time-reversal  imaging  method 
uses  a  double-pass  scanning  scheme  in  a  multipath-rich  environment.  Finally,  we  have  provided 
experimental  results  for  several  key  aspects  of  time  reversal,  such  as  spatio-temporal  focusing  and  an 
extended  virtual  aperture. 


High-Resolution  OTHR  Using  Time  Reversal 


19 


Fig.  16  -  Time  reversal  with  polarization-diverse  antennas  (a)  without  multipath  (free  spaee),  and  (b)  with  multipath  interferenee. 
The  eurves  with  red  eireular  marks  represent  vertieal  polarization  and  the  blue  triangular  marks  represent  horizontal  polarization. 
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